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Abstract: We employ the Monte Carlo particle collision code DPMJET3.04 to determine the multiplicity 
spectra of various secondary particles (in addition to 7r°'s) with 7's as the final decay state, that are 
produced in cosmic-ray (p's and a's) interactions with the interstellar medium. We derive an easy-to- 
use 7-ray production matrix for cosmic rays with energies up to about 10 PeV. This 7-ray production 
matrix is applied to the GeV excess in diffuse Galactic 7-rays observed by EGRET, and we conclude the 
non-7r° decay components are insufficient to explain the GeV excess, although they have contributed a 
different spectrum from the 7r°-decay component. We also test the hypothesis that the TeV-band 7-ray 
emission of the shell-type SNR RX J 17 1 3. 7-3946 observed with HESS is caused by hadronic cosmic 
rays which are accelerated by a cosmic-ray modified shock. By the \ 2 statistics, we find a continuously 
softening spectrum is strongly preferred, in contrast to expectations. A hardening spectrum has about 
1% probability to explain the HESS data, but then only if a hard cutoff at 50-100 TeV is imposed on the 
particle spectrum. 



The 7-Ray Production Matrix 

It is generally believed that diffuse galactic 7-rays 
are a good probe of the production sites and also 
the propagation of accelerated charged particles in 
the Galactic plane (see [3] and references therein). 
Exact knowledge of the 7-ray source spectrum re- 
sulting from hadronic interactions is of prime im- 
portance for a proper physical interpretation for the 
diffuse 7-rays. The production of neutral pions 
and their subsequent decay to 7-rays is thought to 
be the principal mechanism for the 7-ray hadronic 
component in cosmic ray interactions, and a num- 
ber of parameterizations have been developed to 
describe 7-ray production in pp-collision [5, 8]. 
On the other hand, in addition to the neutral pion 
production, the direct 7-ray production in cosmic- 
ray interaction was found important [7], as was the 
production of (E ± , S°), {K ± , K°), and 77 parti- 
cles, all of which eventually decay into 7-rays [8]. 
Also, the helium nuclei in cosmic rays and the in- 
terstellar medium are expected to contribute about 
20-30% to the secondary production in cosmic-ray 
interactions, for which the multiplicity spectra may 
be different. In a hadronic interaction, the multi- 



plicity of non-7r° secondaries with 7-rays as the fi- 
nal decay states is about 50% of the 7r° multiplicity. 
These non-7r° secondaries in hadronic interactions 
might present a different 7-ray spectral profile. 

This work presents a careful study of the 7-ray 
production in cosmic-ray (both p and a) interac- 
tions, by accounting for all decay processes includ- 
ing the direct production. For that purpose we em- 
ploy the event generator DPMJET-III [10] to simu- 
late secondary productions in both p-generated and 
He-generated interactions. We include all relevant 
secondary particles with 7-rays as the final decay 
products. For the composition of the ISM, we as- 
sume 90% protons, 10% helium nuclei, 0.02% car- 
bon, and 0.04% oxygen. Around the energy of 7r° 
production threshold, where DPMJET appears un- 
reliable, we apply a parametric model [5], that in- 
cludes the resonance production for the ir produc- 
tion. We thus derive a 7-ray production matrix for 
cosmic rays with energies up to about 10 PeV that 
can be easily used to interpret the spectra of cosmic 
7-ray sources. 

We consider all the decay channels and their decay 
fractions published by the Particle Data Group to 
account for all secondaries (resonances included) 
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calculated by DPMJET and the parametric method. 
The calculation shows the non-7r° resources in 
hadronic interactions have contributed about 20% 
of the total 7-ray photons, mostly from directly 
produced 7-ray s and decays of T), K L and Kg [7]. 

In the cosmic -ray interactions, we calculate the 7- 
ray spectrum density contributed by decays of un- 
stable secondary particles 



nisM 



J2 I dEN CR (E)cl3a(E) 



(in 



(1) 



where -j^p = d ^- 1 (E 7 , E) is the 7-ray decay 
spectrum from secondary species k. Eq. (1) can be 
re-written into 
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thus reducing this problem to a matrix opera- 
tion with the 7-ray production matrix M.. L j for 
which, each element M.^ = ^2 k d ^ E n shows 
the value of the resultant particle energy spectrum 
j%\E^=E,,E C R=E J , with j being the index for the 
generating cosmic-ray particle (p or a) and i being 
the index indicating the 7-ray energy. The energy 
binnings Ei and Ej are defined with good resolu- 
tions [7]. 

We then use the 7-ray production matrix to analyze 
the observed spectra of diffuse Galactic emission 
and of the shell-type SNR RX J1713.7-3946. 

The GeV Excess 

With the 7-ray production matrix, we calculate the 
diffuse 7-ray spectrum generated by the observed 
cosmic-ray spectrum [9]. Fig. 1 shows the ob- 
served GeV-band 7-ray emission from the inner 
Galaxy [4] in comparison with the contributions 
from 7T° decay as well as bremsstrahlung emission 
describe by a power-law spectrum $b (E): 

$b(E) ~ 1.3 x 10 8 7 • — ^ ^ 

V ; O.leV/cm 3 10 22 cm- 2 

j-, x 2.o-r c 
E \ erg 

100 MeV / cm 2 sec sr 



: decay + bremsstrahlung, l~, 

; w„=0.8, N ISU =3E21 — 

a>.=0,4, N IS „=8E21 - 

a).=0.1 , N ISU =3E22 

- ■ EGRET Data, Hunter et al, 1997 
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Figure 1: Diffuse 7-ray spectrum at GeV range, 
shown in comparison with data for the inner part 
of Galaxy at 315° < I < 345° and \b\ < 5°[4]. 
See text for discussions. 

with power-law spectral index T e = 2.1, the elec- 
tron energy density uj e — 0.1, 0.4, 0.8 eV/cm 3 , 
and gas column density Njsm — 3 • 10 22 , 8 ■ 
10 21 , 3 • 10 21 cm -2 respectively. The generating 
cosmic -rays are assumed with an energy density 
Pe = 0.75 eV/cm 3 . Models based on the lo- 
cally observed cosmic -ray spectra generally pre- 
dict a softer spectrum for the leptonic compo- 
nents, even after accounting for inverse Comp- 
ton emission [4], so we may in fact overestimate 
the GeV-band intensity of the leptonic contribu- 
tion. Nevertheless, it is clearly seen in this figure, 
that in the total intensity an over-shooting around 
Ej ~ 300-600 MeV appears in the modelled 7- 
ray energy distribution, whereas a deficit is present 
above 1 GeV. The observed spectrum of diffuse 
emission is always harder than the model spec- 
trum, and we therefore conclude that an inaccurate 
description of hadronic 7-rays is ruled out as the 
origin of the GeV excess. 

TeV-band Emission from RX J1713.7-3946 

For the TeV-band 7-ray spectrum of the shell-type 
SNR J1713.7-3946 observed by the HESS collab- 
oration [6], we use the 7-ray production matrix to 
test cosmic -ray acceleration models [1, 2], which 
predict a continuous hardening of the cosmic-ray 
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Figure 2: The TeV-band 7-ray spectrum observed 
from RX J1713.7-3946 with HESS [6], shown in 
comparison with the best-fit model of hadronic 7- 
ray production of Eq. (4). 



spectrum up to a high-energy cutoff. We therefore 
parametrize the spectrum of accelerated hadrons as 



N(E) = N Q 
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where 6 is the step function and Eq = 15 TeV is 
a normalization chosen to render variations in the 
power-law index s statistically independent from 
the choice of spectral curvature, a. The cutoff en- 
ergy, E maK , is a free parameter. The normaliza- 
tion iVn is obtained by normalizing both the data 
and the model to the value at 0.97 TeV. By the x 2 
statistics, we obtain the best-fitting values and the 
confidence ranges of the three parameters, E max , 
s, and a, given values as s = 2.13, a = —0.25, 
and E m . dx > 200 TeV, i.e., no cutoff. The best 
fit, shown in Fig. 2, involves a continuous soft- 
ening and is thus not commensurate with expec- 
tations based on acceleration at a cosmic -ray mod- 
ified shock [1,2]. Note the need for data in the en- 
ergy range between 1 GeV and 200 GeV that may 
be provided by GLAST in the near future. Fig. 3 
and Fig. 4 shows the confidence ranges of the pa- 
rameters in Eq. (4), for a vs. s, a vs. E mdx and 
s vs. E mdx respectively, with confidence levels of 
1, 2 and 3 sigma. The contour for E mdx are open 
toward higher energies, indicating that a cut-off is 
not statistically required. The analysis shows that 
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Figure 3: The confidence regions for the spectral 
curvature a and the spectral index s. The shaded 
areas correspond to probabilities 68% (1 sigma), 
95% (2 sigma) and 99.7% (3 sigma). 

a very high E mdx is possible, but requires a nega- 
tive spectral curvature, i.e., < 0, with confidence 
more than 95%, in contrast to the expectation of 
standard cosmic -ray modified shock models. 
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Figure 4: Confidence contours as in Fig. 3. 

Note the step function term <d(E m - dx — E) in 
Eq. (4) can be replaced by an exponential term 
exp(— E/E mdx ). Fig. 5 shows the confidence level 
of E mdx vs. s for the new best-fit parametrization. 
The contour shows a minor deviation from original 
confidence contours in Fig. 4, however, the conclu- 
sion is not altered. 
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Figure 5: Confidence contour as in Fig. 3 but for an 
exponential term exp(— E/E m!ix ) in Eq. (4) instead 
of the step function. 

Conclusions 

We have considered a full picture of the hadronic 
7-rays in cosmic -ray interactions and introduced 
an easy-to-use 7-ray production matrix which can 
be used for arbitrary cosmic -ray spectrum. The 
matrices are available for download at website 
http.V/cherenkov.physics. iastate. edu/gamma-prod. 
We apply the production matrix to calculate the 
7-ray GeV excess and also the TeV-band spectrum 
of SNR RX J1713.7-3946. We conclude that 1) 
the modifications in the GeV-band 7-ray emission 
of hadronic origin are insufficient to explain the 
GeV excess in diffuse galactic 7-rays; 2) a soft 
cut-off at about 100 TeV is statistically required 
in the particle spectrum if the TeV-band spectrum 
of RX J 17 13.7-3946 as observed with HESS is 
caused by cosmic-ray nucleons; 3) no evidence 
for efficient nucleon acceleration to energies near 
the knee in the cosmic-ray spectrum, nor evidence 
of the spectral curvature and hardness predicted 
by standard models of cosmic-ray modified shock 
acceleration. We emphasize the need for GLAST 
data to better constrain the 7-ray spectrum below 
100 GeV. 

In addition to 7-ray production matrix, we have 
also produced the production matrices of all par- 
ticles that are stable on a time scale relevant for 
cosmic -ray propagation. The leptonic production 
by the same approach will be discussed in an inde- 
pendent article. 



Acknowledgments 

Grant support from NASA with award No. NAG5- 
13559 is gratefully acknowledged. 

References 

[1] E. Amato and R Blasi. Non-linear particle 
acceleration at non-relativistic shock waves 
in the presence of self-generated turbulence. 
MNRAS, 371:1251-1258,2006. 

[2] E. G. Berezhko and D. C. Ellison. A Simple 
Model of Nonlinear Diffusive Shock Accel- 
eration. APJ, 526:385-399, 1999. 

[3] A.D. Erlykin, M. Lipski, and A.W. 
Wolfendale. High energy cosmic ray 
spectroscopy. IV. The evidence from di- 
rect observations at lower energies and 
directional anisotropies. Astropart. Phys, 
8:283-292, 1998. 

[4] S. D. Hunter et al. EGRET Observations of 
the Diffuse Gamma-Ray Emission from the 
Galactic Plane. APJ, 481:205-240, 1997. 

[5] T. Kamae et al. Parameterization of 7, e ± 
and Neutrino Spectra Produced by p-p Inter- 
action in Astronomical Environments. APJ, 
647:692-708, 2006. 

[6] F. A. Aharonian et al. The HESS Collabora- 
tion. A detailed spectral and morphological 
study of the gamma-ray supernova remnant 
RX J1713.7-3946 with HESS. AA, 449:223- 
242, 2006. 

[7] C.-Y. Huang, S.-E. Park, M. Pohl, and CD. 
Daniels. Gamma-rays produced in cosmic- 
ray interactions and the TeV-band spec- 
trum of RX J1713.7-3946. Astropart. Phys, 
27:429-439, 2007. 

[8] S. R. Kelner, F. A. Aharonian, and V. V. 
Bugayov. Energy spectra of gamma rays, 
electrons, and neutrinos produced at proton- 
proton interactions in the very high energy 
regime. PRD, 74:034018-034033, 2006. 

[9] M. Mori. The Galactic Diffuse Gamma-Ray 
Spectrum from Cosmic-Ray Proton Interac- 
tions. APJ, 478:225-232, 1997. 
[10] S. Roesler, R. Engel, and J. Ranft. Advanced 
Monte Carlo for Radiation Physics. In Par- 
ticle Transport Simulation and Applications 
(MC 2000), 2000. 



